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On the Application of Some Interpolating 
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J. M. Whittaker considered interpolating equidistant samples by using the cardinal series. In 
this paper we consider the extension provided by the generalized sampling theorem of Kramer for 
nonequidistant samples. We calculate some sampling functions for cases of interest in mathematical 
physics. In particular, this includes the S-matrix contribution due to Regge poles, especially when a 
series expansion other than that of Fourier-Legendre type is needed. 
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1. Introduction 

J. M. Whittaker [1] l considered the interpolation of functions sampled at equidistant points. 
This resulted in his cardinal series and its integral representation for the functions as finite Fourier 
transforms. The result was introduced to Communications Theory by Shannon [2] in the form: if 

fW=J° a e ixt g(x)dx, (1) 

where 

g(x)eL 2 (— a, a), then 

f(t) = ±f(^ Y hi(at - nn) . (2) 

J. M. Whittaker came verv close to touching the question of interpolating functions sampled 
at nonequidistant points when he considered a more general partial fraction series with samples 
at {C,,}, a strictly increasing sequence of positive numbers which is square summable. He noted 
that his cardinal series and its finite Fourier transform representation did not apply but hinted 
that a similar theorem would hold if C t1 =j {)M , the rath zero of the Bessel function Jo{x). So it is 
no surprise to find Bessel functions among the first examples of the following Kramer sampling 
theorem [3]. 

Let / be an interval. Suppose that for each real t 

f(t)=\K(t,x)g(x)dx, (3) 

where g(x) e l 2 (I). Suppose also that for each real t, K(t, x) e L 2 (I) and that there exists a countable 
set E = {t n } such that {K(t„, x)} is a complete orthogonal set onL 2 (J). Then 

f(t) =lim 2 f(tn)S n (t), (4) 

N-kx n\<N 
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/here 



Sn(t) 



K(t, x)K{t n , x)dx 
f \K(t n ,x)\ 2 dx 



(5) 



Kramer showed that these conditions on the kernel K(t, x) are satisfied by the solution of any nth 
order, self-adjoint differential equation and gave as illustrations a first-order equation and the 
Bessel equation. Before Kramer, Weiss [4] had arrived at the same result for n = 2. Kramer illus- 
trated this theorem for the case of K(t, x) =J m (xt) with the sampling function 



Sn(t)=- 



2>Jm,nJ m(t) 



J m\Jm, n) O 



(6) 



Urn, n-t 2 ) J m+\(jm,n) ' 

and Campbell [5] illustrated it for the case of K(t, x) = Pt-i/2 (x) , the Legendre function, to obtain 

2n+ 1 sin (t — n— i) 



S n (t)=- 



ir {t—n—i)(t+n+i) 



(7) 



In the next section we calculate S n (t) for the case of the associated Legendre functions, the 
Gegenbauer functions, a Chebyshev function, the prolate spheroidal functions and a Bessel-like 
function satisfying a fourth-order differential equation. For more details of the calculations one may 
consult [6]. 

2. Further Application of the Kramer Generalized Sampling Theorem 

(i) Associated Legendre functions. Consider the Legendre equation [7, p. 121]: 



d*u 



du 



■ l -^dT*- 2z d z 



r 



»(*+!)- j 



u = 



where z, t, fi are unrestricted. The solution is given by: 

1 /z+lY 12 



u(z)=PHz)- 



r(i-M)\* 



1\ M/2 / 1 z\ 

YJ 2 Fi(-M+1;1-m; 2~2J' |1-2|<2 - 



(8) 



(9) 



To obtain the sampling function S n (t) in (5) we need the integral I Pf(x)P™ {x)dx, where n and 

m are integers and t is unrestricted. We substitute the infinite series expansion for P™(x) [7, p. 167] 
and integrate term by term to obtain 



S n (t)=(-l)> 



2n-\-l {n — m) ! Y(t+m+\) sin ir(t-n) 



77- (n+m) \T(t—m+l) (t — n)(t + n+l) 
(ii) Gegenbauer functions: Consider the differential equation [7, p. 178]: 



(10) 



A 

dz 



il-2 2 )(l-2 2 )"- 1 ' 2 



du 

~dl 



+ \,(l-z 2 ) i -'% =0, k t = t(t + 2v), 



(ID 



satisfied by the Gegenbauer functions C v t (x). For S n (t) in (5) we need 



/> 



(x)Q l (x){l-x 2 y- 1 ' 2 dx, 



t being an unrestricted real parameter. We use [7, p. 159, eq (27), p. 178, and eq (23), p. 140] to 
obtain 

m 



2 m mUz 2 — 1) 2 1 

Cr l Hx) = ( 2ot) ; P&M , m + { ->Q. 



(12) 
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If we combine this result with (10) we obtain 



S„(t)=- 

77 



2 2 "'ml 
(2m) I 



2 tt!(tt + m + ^)r 2 (m + j)rU + 2m-hl) smnjt-n) 
(2m + rc)!(f + n + 2m + l)T(t+ l)ir(t-n) 

1 



(13) 



(iii) Chebyshev functions. These are defined as T>v(x) =Cf(x) and U t (x) =C}(x) but we cannot 
use (13) for either. To treat the latter case we write 

2M-i/2r(* + 2/x)r (,/i+i) l=*z 

C?(x) = - r(2/x)r( ^ + Y) (1 ~* 2) 4P ^-"^W, Re^>0; (14) 

compare (12), and use [7 p. 169, eq (1)] 

£ PfUW*)^ {t _ aH ] +(r+1) * [*(*-er)Py(*?Pfr(*) + (o-+|t)ff (x)Pg-iU) (15) 



to obtain 



P [/(W^Wd-^^/i^f^CjWCUiXi-x 2 ) 1 ^ 



1 



2r(* + 2)T( 



' + "KD" 



(t-n)(f + ra + 2) [r(2)] 2 r(<+l)n! 

+ nPfHh MPn% (x) - tPrJi% (x)Pd9k (x)]U. 

This result plus the orthogonal property 



[x{t-n)P~^ (x)P~\f, (*) 



I 



{U n {x)) 2 (\-X i )^ i dx = TT 



(16) 



yields the required S n (t). 

The result for T t (x) = C?(x), the Chebyshev function of the first kind is not available, since 
the above treatment is limited by the condition Re /jl > in (14). 

(iv) Spheroidal wave functions. Consider the differential equation [8, p. 134, eq (1)]: 



u-*>S-*4 + 



\ + 40(l-z 2 ) 



P 



y=0, 



(17) 



where X, 0, jjl are given real or complex parameters, z is a complex variable and the eigenvalue 
X^(6) =t(t + l) at 6=0. For fji=m, an integer, a solution of (17) which remains bounded at z= 1 
is the prolate spheroidal function, Psf(x, 6). For calculating S n '(t) in (5) we need the value of the 

integral I Psf (*, 6) Psf, (x, 8)dx. First we substitute [8, p. 138, eq (22)], 



Psf(x,d)= J (-l) r a?» r (0)Pr + 2r(*), 



(18) 
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then integrate term by term, using the infinite series expansion for P^ 2r (x) in terms of P%(x) 
[7, p. 167] and the result [8, p. 158, eq (1)], 



/: 



(-l) r a$, r (0) 



(n±m)\ 



(n + i)(n-m)\ 



n — n =2r, 



Ps%(x)P™(x)dx= { 



(19) 



0, 



n — n negative or odd, 



to obtain 



/: 



' P<">,(r ft)Ps/"(x flWr-2Va'" a'" r(f + 2r+ m+ 1) sin w(t- n') 

fs n \x,9)fsl (x,0)dx ^ a «'.'- a '.T(t + 2r-m + lKf + B'+4r+l)ir(t-»^ 



o 



We note that this new result reduces to the special case quoted in [8, p. 147, eq (6)], since 
if we let t — n — n in (20) we arrive at 

J-i r 4l (/i +2r— m)! /i + 2r + i 



and when [8, p. 147, eq (7)] is used this becomes 



/: 



[PsV(x)] 2 dx-- 



1 (tt+/7l)! 

Ai + i (ai — /?i) ! 



2r^m-n', (21) 



(22) 



To evaluate S n (t) in (5) we merely substitute (20) and (22). 

(v) Fourth-Order Differential Equations. The following self-adjoint boundary value problem 
is of interest in signal detection [9] 

[D 4 - (m 2 + m 2 )D 2 + m 2 m 2 ]y(te-*) = t 2 e- 2x y{te~ x ) . 

r (*)=0,y'(* n )=0, rc=l,2, . . . . (23) 

If we take a new variable u — e~ x , the solution H m {t n u) is a Bessel-like function which satisfies 

Kramer's conditions. To obtain S n (t) for (5) we need the integral I uH m (tu)H m (t n u)du. For this 

we consider (23) with t replaced by t n , and eliminate the constant coefficient terms between this 
equation and the original (23) to arrive at 



S n (t) 



2tH' m (t) 



(t 2 -t 2 )H'; n (t n ) 



3. Applications in Physics 

The scattering amplitude of two spinless particles interacting via a spherically symmetric 
potential may be written in the form [10] 



1 



A(k,z)= — ^(2l+l)[S l (k)-l]P l (z), 



(25) 



where k is the wave number, z — cos 0, Pi (x) is the Legendre polynomial, Si is the /th eigenvalue 
of the unitary S-matrix, and 



Si— 1 



= ik J A(k,z)Pi(z)dz. 
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(26) 



When the A(k, z) of (25) is represented by the Watson transform as a contour integral in the 
A-plane, where A represents the complex angular momentum, it is found that the contribution 
to A (A:, z) from the nth Regge pole of the S-matrix at \ = a n with residue ft,, is 

y(a„ + i)/3 B P „(-z)/sin (ira„). (27) 

To obtain the SV ,— 1 corresponding to this term, eq (26) may be used. Hence we need to 
evaluate the integral I l } a„(— z)P/{z) dz, which is the special case of (5) that corresponds to the 
Legendre function P a „(x) = K(a„, x). This integral is given in (7) and was used in [10] to give 

St»> - 1 = Mggn+1) § (28) 

(/-««) (Z+a,.+ l) 

We point out that at times it is necessary to obtain series expansions for A(k, z) in terms 
of other orthogonal functions such as the associated Legendre polynomials, the Gegenbauer 
functions, the Chebyshev polynomials, the prolate spheroidal functions and the Bessel-like func- 
tions. In such cases a Watson-like transform corresponding to each of these orthogonal functions 
may be defined in a similar way. Furthermore, for each of these transforms, the S\ n) — 1 term 
corresponding to the rath Regge pole contribution to A{k, z) may be obtained in a way similar 
to that used for (28). This calculation may be carried out by using the sampling functions corre- 
sponding to the orthogonal functions in the kernel of the transforms. These were derived in (10), 
(13), (16), (20), and (24). 
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